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Cosmology In Terms Of The Deceleration 
Parameter. Part II 


Yu.L. Bolotin, V.A. Cherkaskiy, O.A. Lemets 
D.A. Yerokhin and L.G. Zazunov 

’’All of observational cosmology is the search 
for two numbers: Hq and go-” 

Allan Sandage, 1970 


Abstract 


In the early seventies, Alan Sandage defined cosmology as the search for 
two numbers: Hubble parameter Hq and deceleration parameter go- The first 
of the two basic cosmological parameters (the Hubble parameter) describes 
the linear part of the time dependence of the scale factor. Treating the 
Universe as a dynamical system it is natural to assume that it is non-linear: 
indeed, linearity is nothing more than approximation, while non-linearity 
represents the generic case. It is evident that future models of the Universe 
must take into account different aspects of its evolution. As soon as the 
scale factor is the only dynamical variable, the quantities which determine 
its time dependence must be essentially present in all aspects of the Universe’ 
evolution. Basic characteristics of the cosmological evolution, both static and 
dynamical, can be expressed in terms of the parameters Hq and go- The very 
parameters (and higher time derivatives of the scale factor) enable us to 
construct model-independent kinematics of the cosmological expansion. 

Time dependence of the scale factor reflects main events in history of 
the Universe. Moreover it is the deceleration parameter who dictates the 
expansion rate of the Hubble sphere and determines the dynamics of the ob¬ 
servable galaxy number variation: depending on the sign of the deceleration 
parameter this number either grows (in the case of decelerated expansion), 
or we are going to stay absolutely alone in the cosmos (if the expansion is 
accelerated). 

The intended purpose of the report is reflected in its title — ’’Cosmol¬ 
ogy in terms of the deceleration parameter”. We would like to show that 
practically any aspect of the cosmological evolution is tightly bound to the 
deceleration parameter. 

It is the second part of the report. The first part see here http://arxiv.org/abs/1502.00811 



Contents 


1 Deceleration Parameter in Different Cosmological Models 2 

1.1 SCM. 

1.2 Cosmology with power and hybrid expansion laws. 

1.3 Cosmological models with constant deceleration parameter (Berman 

model). 

1.4 Linearly varying deceleration parameter. 


2 Kinematic Aspects of Dynamical Forms of Dark Energy 

2.1 Quintessence. 

2.2 Chaplygin gas. 

2.3 A'-essence. 

2.4 Polytropic model. 

2.5 Tachyon. 

2.6 Chameleon fields. 

2.7 Phantom fields . 

2.8 Crossing the Phantom Divide. 

2.9 Quintom Models . 

2.10 Interaction in the dark sector . 

2.10.1 Model of Universe with time dependent cosmological 

constant . 

2.10.2 Varying A as a function of the scale factor. 

2.10.3 Interacting DM and DE. 


0 

14 

21 

24 

29 

32 

M 

38 

41 

45 

50 


51 

56 


1 

























Chapter 1 

Deceleration Parameter in 
Different Cosmological Models 


1.1 SCM 


We now consider the evolution of the deceleration parameter (DP) in the 
standard cosmological model (SCM) p]-|5].We recall that in the Big Bang 
model, the substances that filled up the Universe - matter and radiation - 
could only provide decelerating expansion. According to the SCM, presently 
the Universe is dominated by the dark energy (DE) - the component with 
negative pressure. It is that component which leads to the observable accel¬ 
erated expansion. Let us determine the redshift and transition time to the 
accelerating expansion, i.e., find the inflection point in the curve describing 
the time dependence of the scale factor. In the SCM, the second Friedmann 
equation 

a 1 

a = Mf} P 

can be reduced to the form 

— = —Hq [2Dao — D m0 (l + ^) 3 ] (1-1) 

CL Z 

where Dao and h2 m0 are the present-day values of the relative DE density 
in the form of a cosmological constant and of the matter density. Hence, 
we obtain the redshift value z t at which the transition from decelerating to 
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accelerating expansion took place as 


2^A 0 y /3 

o _ / 


( 1 . 2 ) 


For the SCM parameters f2 A o — 0.73 and O m0 ~ 0.27, we have z t ~ 0.745. We 
note that result (11.21) can be obtained using the fact that for two-component 
SCM Universe the DP is 


1 _ 3 fl A o 

2 2 (1 + z) 3 Q m o + D ao 


(1.3) 


The condition q = 0 allows to reproduce (II.2ft . 

Let us find asymptotes of the DP in the SCM. For early Universe (z —>■ 
oo) filled with positive pressure components q(z — > oo) = 1/2, i.e. the 
expansion is decelerated as was expected, while in far future with domination 
of the cosmological constant one observes the accelerated expansion q(z —y 
—1) = —1. The latter result is a trivial consequence of the exponential 
expansion of Universe a(t) oc e Ht in the case of domination of DE in form of 
the cosmological constant. 

Note that the dependence q(t ) can be obtained immediately from the 
definition q = —a/(aH 2 ) making use of the SCM solutions for the scale 
factor 


a{t) = A 1 / 3 sinh 2/,3 (t/t A ); 
iZ A0 o 


As the result one obtains 


Q(t) 


1 

2 


1 — 3 tanh 2 ( — ] 

WJ 


(1.4) 


(1.5) 


The dependence q{t) is presented on Fig ll.ll 

Let us now define the time corresponding to the transition to accelerated 
expansion. Inverting (jl.4|) one obtains 


t(a) 


2 _ 


-D 


!/2 TT-l 


A 0 


H 0 arsli 



1/2 

a 3 / 2 


( 1 . 6 ) 
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Figure 1.1: Time dependence of the DP in the SCM. 
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Transforming from the redshift to the scale factor 


a* = (1 + z*) 


*\—1 


Q 


mO 


20 


AO 


1/3 


one finds 

t* = t(a*) = -O^Q //2 i/ 0 ' 1 sinh _1 (l/2) ~ 5.25 Gyr. (1.7) 

o 

Because of physical importance of the obtained result let us give one more 
(seemingly the simplest) its interpretation. If the quantity aH = a grows 
then a > 0 which corresponds to accelerated expansion of the Universe. 
According to the first Friedmann equation 


aH I o 3 Oao + O. m0 /0.73a 3 + 0.27 

~H o V a ~ V a 


( 1 . 8 ) 


It is easy to show that the function above starts to grow at a* = 0.573, 
which corresponds to z* = 0.745. It is worth noting that the transition 
to accelerates expansion of the Universe (z* = 0.75) occurred considerably 
earlier then the transition to the DE domination epoch (z* = 0.4). 

A natural question arises: what is limiting value of the equation of state 
(EoS) parameter w which still provides accelerated expansion of the Universe 
at present time. As we have seen above, the condition of the accelerated 
expansion reads 

^2(Pi + 3 Pi) < 0. 

i 

In the SCM this condition transforms to 

w de < w de < 0.46. (1.9) 


Of course the substance with such EoS is distinct from the cosmological 
constant and van be realized for example with the help of the scalar fields. 

Let us now use the SCM parameters to estimate the absolute value of the 
cosmological acceleration. Taking derivative of the Hubble law w.r.t. time 
one obtains 

V = (H + H 2 )R. (1.10) 

Time derivative of the Hubble parameter 


da — a 2 




a 


(l.n) 
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Then 


H+H =- = 
a 


6 M 


2 (Pm 2p A ) 3m2 


pi 


pi 


Pa g P m 


= h 2 n 


2 


Finally one gets for the acceleration V an analogue of the Hubble law 


( 1 . 12 ) 


V = HR; H 


H 2 



(1.13) 


For instance, at distance R = 1 Mpc one finds V ~ 10 _11 cm/sec 2 . 

Let us now turn back to the hypothesis about equality to zero of the 
average DP (see Part I, section 2.5), expressed by Lima [(7■ Over the last 
nine years of cosmic microwave background observations, the Wilkinson Mi¬ 
crowave Anisotropy Probe (WMAP) results were consistent with the ACDM 
cosmological model in which the age of the Universe is one Hubble time, and 
the time-averaged value of the DP is consistent with zero. As was noted 
in [7], this curious observation has been put forward as a new coincidence 
problem for the ACDM concordance cosmology, which is in fact a ’greater’ 
coincidence than the near equality of the density parameters of matter and 
the cosmological constant. However recent Planck’s results |8] make the new 
coincidence problem not so important. Let us convince ourself in this. Set¬ 
ting t = t 0 and a — a 0 — 1 in (1 1.6 [) . the H 0 t 0 = 1 condition transforms 
to 

e-VW^DAo _ !) + e 3 ' /n ^(D A0 - 1) + 2D A0 + 2 = 0. (1.14) 

The only positive solution to equation (II. 14ft is D A0 fa 0.737125. In general, 
we can write the time-averaged DP (q) as a function of scale factor: 


(?) + 1 


1 

Th 


2 / D m0 ~ i—i /i i 

3 V?W + fW 


(1.15) 


Figure [L2l shows the 95% confidence levels of D A o from both the Planck data 
(0.686 ^o;o4q) and the WMAP 9-year data (0.721 ± 0.050). The probability 
density function of the posterior distribution is indicated by the colorscale. 
Note that q(to) is close to zero only during a brief period in cosmic time. 
Moreover, q(t 0 ) fa 0 in the WMAP 9-year data, however, q(to) % 0 for the 
Planck data. 
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Figure 1.2: Comparing q as a function of the cosmological scale factor, a(t), 
using the WMAP 9-year and Planck data [9]. 
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1.2 Cosmology with power and hybrid ex¬ 
pansion laws 

described by 
(1.16) 


Let us consider a general class of power-law cosmology 
the scale factor 

/ t x “ 

a{t) = a 0 - 
\u 

where a is a dimensionless positive parameter. The DE in this case 


Q = 


a =I-i. 


(1.17) 


aH 2 a 

The positivity of a leads to q > —1. The scale factor in terms of the DP is 


a(t ) = a 0 


t \ -J+ 1 

to. 


(1.18) 


We observe that q > —1 is the condition for expanding Universe in the power- 
law cosmological model. The expansion history of the Universe in power-law 
cosmology depends on the parameters H 0 and qo 


H(z) = H 0 (l + z) 1+q . 
The comoving distance (1406.2445) 


(1.19) 


r(z = 


n c = 


F(x) = 


VWc 


,F \\n c \ dz' 


Hn 


H(z') 


a ; 


2 H 2 ' 

o^o 


sinh(x), k = — 1 
x, k = 0 

sin(af), k — 1 


with substitution (11.201) can be transformed to 

1 


r\z) = 


\/ia 


=F(in c r i [i-(i+ Z )-»]). 


( 1 . 20 ) 


( 1 . 21 ) 











Note that the Milne model can be viewed as a power-law cosmology with 
ft c = l,a = lorg = 0. For this model therefore 


r(z) = sinh[ln(l + z)\. 


( 1 . 22 ) 


In power-law cosmology the scale factor a(t) and the CMB temperature T(t) 
are related through the relation 


Consequently, 



(1.23) 


(1.24) 


Let us consider now a simple generalization of power-law cosmology, called 
the hybrid expansion law [ 13] 


Cl(t) — <2q 



(1.25) 


where a and (3 are non-negative constants. Further, ao and to as well as in 
the power-law cosmology respectively denote the scale factor and age of the 
Universe today. The cosmographical parameters : Hubble parameter, DP 
and jerk parameter are respectively 


LXb o _ 2 - 

(,fit + at 0 ) 2 

atl(2t 0 — 3 (3t — 3ato) 

((3t + ato ) 3 + 


(1.26) 


In particular cases, one obviously obtains power-law and exponential expan¬ 
sion from (11.25ft choosing a = 0 and /3 = 0 respectively. Accordingly, for 
t —> 0, the cosmological parameters approximate to the following: 

f t\ a TT a 1 32 . 

a —)• ao ( — ) , H —> — , q —> —1 4—, j 1- 1—(1-27) 

V to J t a ol a z 
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Similarly, the exponential term dominates at late times, such that in the 
limit t — > oo we have 

a —)■ a oe^o 1 ), H —> —, q —> — 1, j —y 1. (1-28) 

It may be observed that the parameter a determines the initial kinematics of 
the Universe while the very late time kinematics of the Universe is determined 
by the parameter f3. When a and (3 both are non-zero, Universe evolves 
with variable DP given by (ll.26p and the transition from deceleration to 
acceleration takes place at 

t y/oi — OL 

which puts a in the range 0 < a < 1. 

1.3 Cosmological models with constant de¬ 
celeration parameter (Berman model) 

In 1983, Berman [2] proposed a special law of variation of Hubble parameter 
in FLRW space-time, which yields a constant value of DP. Such a law of 
variation for Hubble’s parameter is not inconsistent with the observations 
and is also approximately valid for slowly time-varying DP models. The law 
provides explicit forms of scale factors governing the FLRW Universe and 
facilitates to describe accelerating as well as decelerating modes of evolution 
of the Universe. Models with constant DP have been extensively studied in 
the literature in different contexts [T51 - 120] . 

Let us consider a law of variation for the Hubble parameter proposed by 
Berman, which yields a constant value of DP, 

H = Da~ n . (1.29) 

Integration of this relation leads to 

a(t) = (nDt + ci)™. (1.30) 

where ci is a constant of integration. The value of DP in this case is 

q = n- 1 (1.31) 
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which is a constant. As usual, sign of q determines the expansion regime. A 
positive sign of q , i.e., n > 1 corresponds to the standard decelerating model 
whereas the negative sign of q, i.e., 0 < n < 1 corresponds to accelerating 
model. The expansion of the Universe at a constant rate corresponds to 
n — 1, that is, q = 0. Therefore in the considered model one can choose the 
values of DP consistent with the observations. 

1.4 Linearly varying deceleration parameter 

We have considered above the simplest case — the time-independent DP. 
Such situation clearly contradicts to our understanding of the Universe’s 
evolution: changes of epochs must inevitably lead to variation of the DP. 
Recall that in the standard ACDM model, the DP is variable and evolves 
from 1/2 to —1. Linear parametrization of the DP represents quite naturally 
the next logical step. In particular, the parametrizations q = q 0 + q^z and 
q — Qo + <7i(l — 0,0 /a) of DP, which are linear in cosmic redshift z and scale 
factor a, have been frequently utilized in the literature [2T1 fI5\ . 

A more general approach is to expand the DP in Taylor series 



(1.32) 


where x is a some cosmological parameter as cosmic scale factor a, cosmic 
redshift z, cosmic time etc. The next step is to choose the parameter to 
decompose in and number of the terms kept. First of all the choice is dictated 
by the adequate description of the observational data, which of course will 
improve when the number of terms in Taylor expansion is increased. However 
increasing the number of terms will lead to growth of the model parameter 
number and make it less efficient. Decreasing of the decomposition terms 
number can be possibly compensated by a lucky choice of the cosmological 
parameter to decompose in. Therefore as the first step in utilization of the 
decomposition (jl.32j) (it is natural to call this procedure the parametrization 
of the DP) one can take the following linear approximation 



(1.33) 


Detailed analysis of the linear approximation III.32)1 was performed in [261127] . 
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Let us start from consideration of the linearly varying DP in terms of 
cosmic redshift z. In this case 

x a 0 

— = — = 1 + z 
xo a 

and 

q = Qo + qiz. (1.34) 

We note that the DP grows monotonically with no limits as we go back to 
earlier times of the Universe. On the other hand, the model is well behaved 
in the future such that the Universe either expands forever (provided that 
qo > qi ~ 1) or ends with a Big Rip in finite future (provided that qo < 
qi — 1). However, the predicted future of the Universe using (jl .34(1 cannot 
be considered so reliable since the observational constraints obtained using 
it with high redshift data are not reliable. 

Unlimited growth of the DP for large z in the parametrization (jl.34j) 
forces us to consider the linearly varying DP in terms of scale factor which 
can be obtained by substituting 

x a 

Xq (Iq ’ 

q = qo + qi (i - —V (1-35) 

V a oJ 

Such parametrization can be written in terms of redshift as follows: 

q = qo + qi 7 ^—- (1-36) 

Treating the Universe as a dynamical system it is useful, along with the 
parametrizations fll.34p - fll.36p . to consider the parametrization of the DP 
directly in terms of cosmic time t. Such parametrization can be obtained by 
substitution 

x t 
Xo to 

in QEE): 

? = 9o + ?ifl-^- (1.37) 

In contrast to the two previous models, this model cannot be used for ob¬ 
servational analysis directly. One should first solve it for the scale factor 
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explicitly, and then obtain the time red-shift relation. This procedure can be 
facilitated if one first solve the equation with separating variables 


dH 

dt 


Substituting (11.371) . one finds 


- H 2 (l + q ). 


H(t) — 2—[2(1 + qo + qi)to — cjit] L . 


Using the obtained expression for H(t) to integrate the equation 

dz , 

- = -(l +z)H(t) 

we obtain: 


q = qo + qi 



2('\ 4- On 4- «i ) 


(1.38) 


(1.39) 


(1.40) 


(1.41) 


There are many ways to parametrize the DP with different forms of EoS 
for components filling the Universe. Let us consider for instance a single¬ 
component flat Universe with a fluid described by an EoS parameter ex¬ 
pressed as a first order Taylor expansion in cosmic time [25], that is: 


W = W 0 + M7i(l - t), 


where Wq and W\ are real constants and t is the normalized time. This 
parametrization in the single-component flat (0 = 1) case gives for the DP 

q 3 _ ^ 3 

q= 2 + 2 = 2 + 2 ^ Wo + Wl ^ “*)]■ 

i 

Such an approach enables us to use the cosmographic (model-free) constraints 
to test the EoS for different components. 
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Chapter 2 

Kinematic Aspects of 
Dynamical Forms of Dark 
Energy 

2.1 Quintessence 

The cosmological constant is just one of possible realizations of the hypothet¬ 
ical substance — the DE, introduced to explain the accelerated expansion of 
the Universe. As we have seen above, such a substance must have the EoS 
p = wp with parameter w which satisfies the condition w < — 1/3 (in absence 
of other components). Unfortunately the nature of the DE is mysterious for 
us and that is why there is huge number of hypotheses and candidates for 
this role of the fundamental component of energy budget of the Universe. 
We mentioned many times about the fast progress of observational cosmol¬ 
ogy during the last century. However we still cannot answer the question of 
time evolution of the parameter w. If this parameter does depend on time 
then we have to look for an alternative to the cosmological constant. For 
a short time a great number of the alternative possibilities (w.r.t w = —1) 
was studied. Scalar fields which formed the post-inflationary Universe are 
one of the main candidates for role of the DE. A version of scalar field p 
with appropriately chosen potential V(p) is the most popular one. In such 
models, in distinction to the cosmological constant, the scalar field is a dy¬ 
namical variable, and density of the DE depends on time. The models can 
be distinguished by the choice of Lagrangian for the scalar field. Let us start 
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from probably the simplest model of the DE of such type called quintessence. 

The quintessence means the scalar held p> in the potential V ( ip ) minimally 
coupled to gravity, i.e. affected solely by the curvature of space-time, and 
limited by canonical form of kinetic energy. Action for such a held reads 


S 


J d A x\/—gL = J d A x\J—g 


1 nv 9 V d <P 

2 ^ dx^ dx v 


V(<p) 


( 2 . 1 ) 


where g = det g tlv . Equation of motion for the scalar held is obtained by 
variation of the action w.r.t. the held 


1 

7=9 


9 v 




dV(ip) 

dip 


( 2 . 2 ) 


In hat Friedmannian Universe, i.e. for the FLRW metric, one obtains for the 
homogeneous held <p(t) 


ip + 3 Hip + V'(<p) = 0. 


(2.3) 


For the case of homogeneous held ip(t) in locally Lorentzain frame, where the 
metric g llv can be replaced by the Minkowski one, one obtains the following 
expression for density and pressure of the scalar held 

P v = Too = - ip 2j rV(<p)] (2.4) 

P<p = Ta = ^0 2 -U(v7). 

The Friedmann equations for hat Universe filled with the scalar held, and 
the Klein-Gordon equation take on the form 


H 2 = - 


H — — 


hp 2 + V(ip) 

t. 

2 ’ 


dV 

p> + 3 Hip -|——— = 0. 
dip 


(2.5) 


Using (12.4p one obtains the EoS for the scalar held 













The EoS parameter w v for the scalar field varies in the interval 


— 1 < w v < 1 . 


(2.7) 


The EoS for the scalar field can be conveniently presented as 

, , x — 1 


W v {x) = 


X + 1 


^ 2 /2 


( 2 . 8 ) 


The function w(x) monotonically grows from the minimal value w m i n = — 1 
at x = 0 to maximal asymptotic value w max at x —* oo, which corresponds 
to V = 0. In the slow-roll limit x <C 1 (<p 2 <C V(tp)) the scalar field behaves 
as the cosmological constant w v — —1. It is easy to see that = const 

in this case. In the other limit x 1 (<p 2 V(<p)) (rigid matter) w v = 1. 
In this case the energy density of the scalar field evolves as p v oc a -6 . The 
intermediate situation x ~ 1, p ~ 0 corresponds to the case of non-relativistic 
matter. 

The DP for the Universe filled with the scalar field is 


= ±(±\_ i _ <?-v 

% dt \H) |0 2 + U‘ 


(2.9) 


Let us consider now a new formalism for the analysis of scalar fields in flat 
isotropic and homogeneous cosmological models, proposed in [29]. Using the 
formalism, the basic evolution equation of the models can be reduced to a 
first order non-linear differential equation. Approximate solutions of this 
equation can be constructed in the limiting cases of the scalar field kinetic 
energy and potential energy dominance, respectively, as well as in the inter¬ 
mediate regime. In the following we will restrict our study to expansionary 
cosmological models, which satisfy the condition that the scale factor is a 
monotonically increasing function of time. For expanding cosmological mod¬ 
els the condition H > 0 is always satisfied. 

From Eqs. (I2.5jl we obtain the Riccati type equation satisfied by H, of 
the form 

H = V — 3H 2 . (2.10) 

By substituting the Hubble function into the Klein-Gordon equation, we 
obtain the basic equation describing the scalar field evolution as 


/-lip 2 dV 

v 3 + v / 3\/ y + U (<p) 0 + — = o. 


( 2 . 11 ) 
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In order to deduce a basic equation describing the dynamics of the scalar 
fields in the flat FLRW Universe, we define a new function G(ip) by 



so that Eq. (I2.11j) becomes 



0 . 


( 2 . 12 ) 


As a function of time G satisfies the equation 


Jr* 

—- = — \fzvjjp) sinh G 



3 1 dV 

2 2U (ip) dip 


coth G 


(2.13) 


Let us now turn directly to the question of our interest — DP as a function 
of the scalar held. Using (12. 12)1 one obtains 



(2.14) 


The latter result evidently coincides with the Eq. ( 12.9j) . 

The considered approach enables us to obtain a number of analytic ex¬ 
pressions for DP in a wide class of scalar held potentials. Let us briehy cite 
below several examples (see [29] for details). 

The exponential potential scalar field. If V'/V is a constant, that is, 
V'/V = Vboo = const , the scalar field self-interaction potential is of the 
exponential form, 



(2.15) 


where Vo is an arbitrary constant of integration, then Eq. (12.12)1 takes the 
form 



(2.16) 


For the case a 0 ^ ±1 Eq. (12. 16)1 gives immediately 



G — «o hi | sinh G + a 0 cosh G \ 



(2.17) 
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where <po is an arbitrary constant of integration. The time dependence of the 
physical parameters can be obtained from Eq. (I2.13P as 



/ 


eX /37i«o^ ( sinh G + a 0 cosh G)' 


dG 


(2.18) 


With the use of Eq. (I2.17p . we obtain the following integral representation 
for the time t, 



/ 


e 1 “o (sinh G + a 0 cosh G) a ° 1 dG, (2.19) 


where to is an arbitrary constant of integration. Thus, Eqs. fl2.14p and f!2.19p 
give a parametric representation of the time evolution of the DP, with G taken 
as a parameter. It is represented, for different values of «o, in Fig. 12.11 For 
the considered range of the parameter a 0 , the Universe starts its evolution 
from a decelerating phase, with q > 0, but after a finite interval it enters in 
an accelerated era, with q < 0. In the limit of large times the Universe enters 
a de Sitter type accelerated phase, with q = — 1. 

Generalized hyperbolic cosine type scalar field potentials. As a 

second example of an exact integrability of the evolution equation, given by 
Eq. (I2.12p . we consider the case in which the scalar field potential can be 
represented as a function of G in the form 



( 2 . 20 ) 


where aq is an arbitrary constant. With this choice, the evolution equation 
takes the simple form 



( 2 . 21 ) 


with the general solution given by 



( 2 . 22 ) 
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Figure 2.1: Plot of the time variation of the DP of the Universe filled with 
an exponential potential scalar held for different values of ao: ao = 1.5 
(solid curve), ao = 2.5 (dotted curve), ao = 3.5 (short dashed curve), 
ao = 4.5 (dashed curve), and ao = 5.5 (long dashed curve), respectively. 
The arbitrary integration constants ipo and U 0 have been normalized so that 
exp (-’■y/3/2a 0 tp 0 'j = y/W Q . 
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Figure 2.2: Time variation of the DP of the Universe filled with a scalar field 
with a generalized hyperbolic cosine self-interaction potential for different 
values of ay ay = 0.1 (solid curve), ay = 0.15 (dotted curve), ay = 0.20 
(short dashed curve), ay = 0.25 (dashed curve), and op = 0.30 (long dashed 
curve), respectively. 


where ipo is an arbitrary constant of integration. The time dependence of the 
DP can be obtained in a parametric form (taking ip as the parameter) as 


t-t 0 = 



dip 


cosh 1+ “i 


(1 + ay) (ip - ip 0 ) 


sinh 


q = 3 tanh 2 \/^ (1 + ay) {ip - ip Q ) 


1 . 


(1 + ay) (ip - ip 0 ) 

(2.23) ' 

(2.24) 


The time variation of the DP is represented for different values of an, in 
Fig. 12.21 In all of the considered models the Universe shows an expansionary, 
accelerated behavior, starting with an initial value q = — 1 of the DP. The 
scalar held potential is increasing in time, leading, in the long time limit, to 
accelerated expansions, with q > — 1. 


20 














Power-law type scalar field potential. A simple solution of the grav¬ 
itational field equations for a power-law type scalar field potential can be 
obtained by assuming for the function G the following form 


G 


arccoth 



a 2 = constant. 


(2.25) 


With this choice of G, the Eq. (12.12p immediately provides the scalar field 
potential given by 


V(v) 


V„ I — 

OL 2 


-2(a2 + l) 



(2.26) 


where Vo is an arbitrary constant of integration. The time dependence of the 
scalar field is given by a simple power law, 


<P(t) 

a 2 


\/2Vo (og + 2) 

OL 2 


1 

“ 2+ 2 l 

(t-to)“ 2+2 • 


(2.27) 


The DP is given by 


q = 2 



V2 V 0 (ot 2 + 2) 
ol 2 


2 

a 2+2 

{t - to) 


2 

a 2 + 2 — 


(2.28) 


2.2 Chaplygin gas 


Multiple attempts to explain the observed accelerated expansion of the Uni¬ 
verse in frames of GR stimulated consideration of a new type of matter as 
candidate for DE. The Chaplygin Gas (CG) realizes one of such possibilities. 
The EoS for CG is 

A 

p= -. 2.29 

P 


Desire to improve the description of the observed dynamics of the Universe 
led to multiple generalizations of the EoS (12.291) . 

The EoS for generalized Chaplygin gas (GCG) [31j is given by 


P- ~ 



(2.30) 


21 





















where 0 < a < 1. The additional parameter a widens possibilities of the 
model. Further modification of the EoS (I2.29p includes a positive term lin¬ 
early proportional to the density 


p = Bp — 


A 

P a 


(2.31) 


Here A, B, a are positive constants and 0 < a < 1. Such a model is called the 
modified Chaplygin gas (MCG). The cosmological models based on the EoS 
(12.29j) and its numerous generalization of the type (I2.30p and (I2.3ip share an 
attractive feature to describe substances with essentially different physical 
properties in frames of single EoS. So for instance the MCG model (12.311) 
for low densities describes a substance generating negative pressure while at 
high density and B = 1/3 this equation describes a radiation dominated era. 

Using the EoS (I2.3ip . we can integrate the conservation equation to obtain 


Pmcg(o) — Po 


Ac 


l-A s 

a 3(l-KB)(l+a) 


1+a 


(2.32) 


where 


Af,= 


A 


1 + B p“ +1 ’ 

B/0 and po is an integration constant. In the considered model of Universe 
one have to introduce only the baryon component f \ (there is no necessity 
for introduction of the DE), as the Chaplygin gas represents an unified model 
of both the dark components. Using 


a 


Q = 


tot 


r\ 'y ^ Wi^li, 


2 2 

l 

for the case of spatially flat Universe one finds 


(2.33) 


q(z) = 


a 3 f2 b0 + VLmcg{°)\\ + 3wMcc( a )\ 


2 [a 3 Q b0 + n MCG (a)] 

^ Mcc ( a ) = ^mcgo[As + (1 — As)a 3 ( 1 +- B )( 1 + Q )] i+a 


(2.34) 

(2.35) 


The EoS parameter w for the MCG can be found using its definition (I2.3ip 

P 


wmcg = - = B - p (q+1) . 
P 


(2.36) 
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(2.37) 


Substitution of the equation (I2.32p into (j2.36j) gives 

w M cc(a) -B- — - (1 _ Ai)a _^ B ^ • 

Let us now express parameters of the CG with the EoS (\2.29\i in terms of 
the cosmographic parameters. As we have seen above for small values of a(t) 
density of the Chaplygin gas reduces to the dust density p oc a -3 , and for 
large a one gets the de Sitter Universe: p = const, p = —p. In between these 
two regimes one can use the approximation 

TD 

p=y /A+^=a- 6 . (2.38) 

y/2A 

Thus \fA plays the role of a cosmological constant. We insert this to the 
Friedmann equation 

a 2 + k = Irpa 2 + ^Aa 2 (2.39) 

o o 

with A and follow the procedure of eliminating the constants by differentia¬ 
tion. This leads to an approximate constraint 


s + 5(q + j) + qj = 0. 


(2.40) 


Analogous procedure for the generalized Chaplygin gas (I2.30p leads to the 
equation 30' 

s + (3a + 2)(q + j) + qj = 0. (2.41) 


in agreement with 

s + 2(q + j) +qj = 0. (2.42) 

and (I2.40p when a = 0 and a = 1 respectively. If we want to exclude the 
parameter a from the latter equation we must take one more derivative of 
the Friedmann equation and introduce an additional cosmological parameter 


1 d 5 a /1 da\ 5 
a dt 5 \a dt J 


As a result one obtains 


— 2qs — 2jq 2 — lq — 2sj — 3sq 2 —j 2 q — lj + s 2 — 3q 2 j—qsj+j 3 — 2j 2 q 2 = 0. (2.43) 

This constraint is again approximate and is valid only in the regime where 
the higher order terms in (I2.38p can be dropped. 
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2.3 iv-essence 


As we have seen above, in the quintessence model the required dynamical 
behavior is attained due to appropriate choice of the scalar field potential. 
Another DE model called the /(-essence [321435] is realized by the scalar held 
with modified kinetic term. 

Let us introduce the quantity 



dp dp 
dx M dx v 


and consider the scalar held action of the form 


S = j d 4 xy/^gF(p,X), 


(2.44) 


where F is generally speaking arbitrary function of the variables p and X. 
Let us hrst consider the p -independent case |36j 


S 


f d l xX=gF(x). 


(2.45) 


Using the energy momentum tensor for a perfect fluid we obtain that the 
fc-essence energy density p and pressure p are 


p = 2XF x — F, p = F. 


(2.46) 


We will assume that the energy density is positive so that 2 XF x — F > 0. 
The EoS is 


p F 

p 2 XF x — F 


(2.47) 


One can use equation (I2.47jl to express the condition w > — 1 (w < —1) as 
a condition on the function F(X). We need to consider the two possibilities 
F > 0 and F < 0 separately. In the hrst case, demanding the energy density 
be positive immediately implies w > 0. For F < 0, a positive energy density 
means that also 2XF x /F < 1 so 


-1 


w = 


1 - 2 XF x /F 


> - 1 , 


(2.48) 
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when F x > 0. Therefore [37 


F x > 0 => w > —1; 
F x < 0 w < —1; 
=>- w > 0. 


F < 0, 
F < 0, 
F > 0, 


(2.49) 


The latter condition excludes the possibility of the accelerated expansion. 

Next, we consider the case F((p, X). The non-canonical scalar held La- 
grangian density is given by Fang et al. [38], 

F(<p,X) = r(X)-V{<p), (2.50) 

where V{<p) is a self-interacting potential for the scalar held <p, T(X) is 
an arbitrary function of X, which is dehned as X = <p 2 /2. Following the 
article [39] let us consider a Lagrangian density of the following form 

F(f, X ) = A' 2 -VM, X - iy. (2.51) 


In this case, the energy density and pressure associated with this Lagrangian 
density read 


/v = + V(ip)-, 

(2.52) 

Pv = 

(2.53) 

and Friedman equations with conservation equations take 

on the form 

3 H 2 = p m + —y? 4 + V(<^); 

(2.54) 

2H + 3H 2 = ~ l p A + V((p); 

(2.55) 

p v + 3 H(p tp Fp^) = 0; 

(2.56) 

Pm T 3 H p m 0. 

(2.57) 


Among this four equations (equations (I2.54I) - (I2.57I) ). only three are indepen¬ 
dent equations with four unknown parameters H, p m , ip and V{ip) and for 
this reason we still have freedom to choose one parameter to close the above 
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system of equations. As for this parameter we can choose parameter of DE 
EoS. Parametrization of EoS parameter as, 

w v (z) = — = w 0 + wxf(z), (2.58) 

Pip 

where Wq, w± are real numbers and f(z) is a function of redshift z. 

The most popular parametrization Chevallier-Polarski-Linder (CPL) parametriza¬ 
tion is one of the most popular one and is given by 


Wip(z) = w 0 + tUi(l - a) = w 0 + W\ ’ 


(2.59) 


where z — 1/a —lis the redshift, w 0 represents the current value of w ip (z) and 
the second term represents variation of the EoS parameter w.r.t. redshift. 
Note that this parametrization has the advantage of giving finite w v in the 
entire range, 0 < z < oo. 

The DP for this model takes the following form 


q(z ) 


i3 r wp+ Wl (^) 

2 2 [i+K(i+^)( 3 - 3 °i)e(^) 


(2.60) 


where 

PmO ^mO 
K= - = n—' 

pcpO ^ 

If q < 0 the model accelerates while q > 0 indicates deceleration of the 
universe. From the figure (12.3h ). we have seen that q decreases from positive 
to negative value by suitable choices of model parameters [39J. 

However, the model presented here is restricted because oj v (z) diverges 
when z —> — 1 , i.e. this model cannot predict about the future evolution. 
So, the model can nicely describe the evolution history of the universe in the 
past and near future up to z > — 1 but can not predict about the evolution 
beyond that limit. 

Recently, Jassal et al. [4Qj extended the above parametrization to a more 
general case: 

( z ) — u 0 + uji + . (2-61) 

In this model the DP is expressed as 


q(z) 


1 

2 


3 

2 


w o + w iJTTW 

r 3™ 1Z 2 \ 

1 + k(1 + z)~ 3w °e VTT T^) 


(2.62) 
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Figure (j2.3b ) shows the plot of q(z) as a function of z. This plot clearly 
shows the transition of q from the decelerating to the accelerating regime at 
z ~ 0.8. In graph, the resulting cosmological scenarios are in good agreement 
with observations. For the JBP model also, co ip (z) diverges as z —» — 1 and 
thus future evolution can not be predicted. 






Figure 2.3: a) Top left: q(z) for CPL parametrization. b) Top right: 
q(z ) for extended CPL parametrization. c) Bottom left: q(z ) for Barboza- 
Alcaniz parametrization. d) Bottom right: q(z ) for Generalized Chaplygin 
Gas parametrization. For all figures n = Ll m0 /Q v0 = 0.27/0.73, w o = — 1 
and w i = 0.2 

The next parametrization was proposed by Barboza et al. [3Tj, which has 
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the following functional form 


w<p{z) = W 0 + Wi 


z(l + z) 

1 + z 2 ’ 


where w v (z = 0) = w 0 (the present value of the EoS parameter), 


w i = 


dw^ 
dz 


z =0 


(2.63) 


(which measures the variation of the EoS parameter with z), w v (z = oo) = 
wo + wi and the EoS parameter reduces to w v (z) = wo + w\z at the low 
redshift (z <C 1). It is remarkable that the BA parametrization does not 
diverge like CPL model when z —> —1 and it remains bounded within the 
entire range z G [—1, oo). 

In this model the DP q(z) can be expressed as 


q(z) 


_ r i z(i+z) 

1 3 Wo+Wl^+72 2 

2 2 1 + r( 1 + ^)- 3,w o(l + x; 2 ) - ^ 1 


(2.64) 


Figure (12.3b ) shows the evolution of the DP with redshift z. We have 
observed from figure (12.3b ) that the evolution of the universe is in accelerating 
phase (q < 0) at present epoch. 

As we have seen in the previous section, the generalized Chaplygin gas 
behaves like dark matter (DM) in the past and it behaves like cosmological 
constant at present. Assuming that the Universe contains both the DM and 
the GCG, we also take into account the possibility that the fc-essence ip plays 
the role of GCG to explore the late time cosmic scenarios. 

By inserting equation (I2.30K into the energy conservation equation (12.56ft , 
one finds that the density of the scalar held <p evolves as 


p v (z) = [A + j B(1 + ^) 3 ( 1+q )]( 1 +“ ) , 


(2.65) 


where B is an integration constant. Equation (I2.65p can be re-written in the 
following form 

pM = Pv 0 |A + (! - 40(1 + z) 3 ( 1+ “)] A), (2.66) 


where, for simplicity, we have defined 


A 


S 


A 

A + B 
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and p v o = (A + B)( 1 + a '> is the present value of the energy density of the 
GCG. To ensure the finite and positive value of p v we need — 1 < a < 1 and 
0 < A a < 1. 


In this case, the Hubble parameter is given by 


H 2 = Hi 


hlmo(l + ^) 3 + (A s + (1 — A s )(l + <1+ “ ) . (2.67) 


The corresponding expression for the EoS parameter is given by 

A. 


W<p(z) = 


H, + (l-H s )(l + ^) 3 ( 1 +“)' 


( 2 . 68 ) 


Like earlier mentioned three models, the EoS parameter of the GCG also de¬ 
pends on two independent model parameters (A s and a) along with redshift 
z. At present epoch, the above EoS parameter becomes, u v (z = 0) = — A s . 
In this case, it is very interesting to notice that the GCG will behave like 
pure cosmological constant when we put A s — 1. 


The DP q can be written as 


_ As _ 

A a + (1-A s )(l+2:) 3 ( 1 + Q ) 

_ k(1+z) 3 _ 

[A s + (l-A s )(l+z) 3 ( 1 +“)] (!+“) 

where 

PmO hlmO 

K = - = —-. 

pipO ^(^0 

Figure (12.3d ) shows the evolution of q(z) vs. redshift z for A s = 0.9 and 
a = —0.5. In fact, at low redshift, the transition of q(z) from decelerating to 
accelerating regime depends upon the choice of the parameters A s and a. 


(2.69) 



2.4 Polytropic model 

Following Chavanis J42l 144] we construct models of Universe with a EoS hav¬ 
ing a linear component and a polytropic component, 

p = ap + kp 1+1/n . (2.70) 
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Concerning the linear EoS, we assume — 1 < a < 1. As for the polytropic 
component, we consider a wide range of parameters. 

With positive index n > 0, the polytropic component dominates the 
linear component when the density is high. These models describe the early 
Universe. Conversely, when n < 0, the polytropic component dominates the 
linear component when the density is low. If 1 + a + kp 1 t n > 0, the EoS 
parameter is w > —1, which corresponds to the quintessential phase. The 
opposite case is 1 + a + kp 1 ^ n < 0, leading to a phantom Universe phase. 
Below we consider behavior of the DP for different values of the parameters 
a, n, k. 

We start consideration from the case 1 + a + kp l t n > 0, n > 0, when the 
equation 

p + 3Hp(l + a + kp 1 ! n ) = 0 
has the following solution 

P* 

P= [(a/a*) 3 P + “)/ n =F 1]"’ 


(2.71) 

(2.72) 


where p* = [(a + l)/|/c|] ?l and a* is a constant of integration. The upper 
sign corresponds to k > 0 (repulsive self-interaction) and the lower sign 
corresponds to k < 0 (attractive self-interaction). 

We can rewrite the EoS (I2.70P as p — w(t)p with 


w{t) 


a ± (a + 1) 



(2.73) 


For a flat Universe 

and consequently 


Q(t) 


1 + 3ct 
2 


± ^(a + 1) 



(2.74) 


(2.75) 


We define a critical density p c and a critical scale factor ac, corresponding 
to a possible inflexion point q = 0 in the curve a(t) 

n/[3(l+a)] 

( 2 . 



pc = 


T 


1 + 3a 


_ 3(1 +a) 


Clr 


T 


1 + 3a 
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Let us now analyze behavior of the DP for positive and negative values of 
the parameter k. We first assume k > 0. In this case the density is defined 
only for a > a*.When a —* a*, q —> oo. When a —* oo, q —> (1 + 3a)/2. For 
a > —1/3, the Universe is always decelerating (g > 0). For a < —1/3, the 
Universe is decelerating when a* < a < a c and accelerating when a > a c . 

We now assume k < 0. In this case the density is defined for all a. When 
a —> 0, q —> —1; when a —> oo, q —» (1 + 3a)/2. For a < —1/3, the Universe 
is always accelerating (q < 0). For a > —1/3, the Universe is accelerating 
when a < a c and decelerating when a > a c . 

Let us briefly analyze the case 1 + a + kp 1 ^ > 0, n < 0. Dependence 
of the density on the scale factor is still described by the formula (I2.71H . 
however now for k > 0 the density is defined only for a < a*. 

In this case again the DP in a flat Universe and inflexion point are defined 
by Eqs. (12. 75[) . fj2.76|) . We first assume k > 0. When a — » 0, q — » (1 + 3a)/2, 
when a —» a*, q —» +oo. For a > —1/3, the Universe is always decelerating. 
For a < —1/3, the Universe is accelerating when a < a c and decelerating 
when a c < a < a*. 

We now assume k < 0. When a —> 0, q —> (1 + 3a)/2, when a —» oo, 
q —» —1. For a > —1/3, the Universe is always decelerating. For a < —1/3, 
the Universe is accelerating when a < a c and decelerating when a > a c . 

As a conclusion of the section we analyze the DP in the phantom Universe 
1 + a + kp x ! n < 0, w(t) < —1. In this case the Friedmann equation (12.711) 
can be integrated into 


P = 


P * 


[1 - (a/a*) 3 ( 1+a )/"] r 


(2.77) 


where, as before, p* = [(a + l)/|/c|] n and a* is a constant of integration. For 
n > 0, the density is defined only when a < a*. 

The DP is defined by Eqs. (I2.73p . (I2.74|) 


q(t) 


1 + 3a 
2 


2( a + - 1 -) 



(2.78) 


A phantom Universe is always accelerating since q < — 1 < 0. For n > 0, 
q — » — 1 when a —> 0 and q —> — oo when a —>• a*. For n < 0, q —» — oo when 
a —> a* and q —> — 1 when a —>■ +oo. 
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2.5 Tachyon 

The tachyon of string theory may be described [45j by effective field the¬ 
ory describing some sort of tachyon condensate which in flat space has a 
Lagrangian density 


L = -y(T)Vl + rWT. (2.79) 


where T is the tachyon field, V(T ) is the tachyon potential and r] lw = 
diag (—, +, +, +) is the metric of Minkowski spacetime. The tachyon po¬ 
tential V(T ) has a positive maximum at the origin and has a minimum at 
T = Tq where the potential vanishes. In Minkowski spacetime the rolling 
down towards its minimum value is described by a spatially homogeneous 
but time-dependent solution obtained from the Lagrangian density 


L = —V V1 — T 2 . 


(2.80) 


During the rolling the Hamiltonian density 


Ph 


V(T) 

VI - r 2 


has a constant value E. Thus 


T = 



(2.81) 


(2.82) 


As T increases V (T) decreases and T increases to attain its maximum value 
of 1 in infinite time as T tends to infinity. Note that as explained in [46] the 
tachyon field behaves like a fluid of positive energy density 



V(T) 

(2.83) 


Vl~T 2 

and negative pressure 


p= -V(T)Vl-f 2 . 

(2.84) 

Thus 


w = - = -(1 -T 2 ) 

P 

(2.85) 
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and therefore, — 1 < w < 0. Note that both the Weak Energy Condition, 
p > 0 and Dominant Energy Condition, p > \p\ hold. However because 


p + 3p 


4m (i - ia . 

y/l - f 2 V 2 / 


( 2 . 86 ) 


the Strong Energy Condition fails to hold for small \T\ but does hold for 
large \T\. 

The discussion above has neglected the gravitational field generated by 
the tachyon condensate. To take it into account we use the Lagrangian 
density 

L = V=S (jjg * V(TW 1 + rtv) ■ (2.87) 

and (12.84(1 for the density and pressure 


167 tG 

In this case the expressions 


remain valid and thus the Friedman equations are 

a 2 k 8 t tG V(T) 


a 2 a 2 3 

a 8 t tG V(T) 


a 


Vl ~ T 2 


ji2 

1 - -T 2 
2 


( 2 . 88 ) 

(2.89) 


It follows from (12.821) that as T increases V(T) increases too, but T could 
decrease since E decreases. In any case T remains positive and so T increases 
monotonically to its maximum value of 1. It follows from (12.881) that if k < 0, 
then a will always be positive. This is because the Weak Energy Condition 
holds, p > 0. From the Raychaudhuri equation (I2.89jl it follows that if 
|T| < 2/3 the scale factor initially accelerates (a > 0), but then, when T 
exceeds y/2/3, the acceleration will stop and deceleration will start. If the 
Universe is flat (k = 0), then a(t) —> constant. In the hyperbolic Universe 
(k — —1) the scale factor increases linearly with time (a —» t). In both cases 
the final state of the Universe is flat, the case k — — 1 being the Milne model. 
In the spherical case (k = 1) the Universe will re-collapse. The possibility of 
cosmic acceleration arises from the positive potential V(T). 


2.6 Chameleon fields 

In all the above considered dynamical models of the DE (quintessence, k- 
essence, phantom fields) the scalar field had only the self-interaction, de¬ 
scribed by the potential V((p). Absence of interaction with other components 
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in the Universe, firstly, looks unnatural, and secondly, it limits possibilities 
of the model. However inclusion of the interaction (the procedure is rather 
simple from theoretical point of view) immediately faces the following fun¬ 
damental difficulty. 

In a local part of Universe (say the Solar System) available for precision 
measurements the theory is agreed with the experiment (observations) by 
introduction of four fundamental forces: the strong, the weak, the electro¬ 
magnetic and the gravitational ones. Introduction of a new type of interac¬ 
tion automatically leads to appearance of the ’’fifth force” which was not yet 
observed. The latter fact imposes strict constraints on the interaction of the 
scalar field with matter: either the interaction must be considerably weaker 
than the gravitational one, or quanta of the field must be too heavy, i.e. the 
interaction must be short-range. 

A natural question arises: is it possible to build a model of DE in form 
of the scalar held interacting with matter which would not contradict the 
equivalence principle well verified on scales of the Solar system? Recall that 
the DE denotes any substance capable to provide the accelerated expansion 
of Universe. Positive answer to the above posed question is given by the 
so-called chameleon models of scalar held [J7H52]. 

The chameleon scalar helds are scalar helds coupled to matter with inten¬ 
sity of order of the gravitational forces (or even higher) and mass depending 
on density of the environment. On the cosmological scales where density is 
negligible this helds are ultralight. On the contrary, for example in vicinity 
of the Earth where density is essentially higher the helds acquire sufficiently 
great mass. In other words, properties of the helds including its value vary 
along with changes in density of the environment. That is why they are 
called the chameleon helds. 

Note that introduction of forces depending on density presents a common 
tendency in physics. Let us cite the following example. In the sixties in or¬ 
der to calculate nuclear characteristics by the Hartree-Fock method different 
types of effective nucleon-nucleon interactions were investigated. It turned 
out that no one among the considered potentials led to consistent quantitative 
description of the nucleus. The problem was solved only after introduction 
of the effective interaction between the nucleons, which depended on density. 

Action for the chameleon held represents sum of Hilbert-Einstcin action, 
scalar held action and a term describing the interaction of the scalar held 
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with matter 


S 


^ - V(<p) + f{p)L r 


\J—gd A x. 


(2.90) 


Unlike the usual Einstein-Hilbert action, the matter Lagrangian L m is mod¬ 
ified as f(p)L m , where f(p) is an analytic f un ction of p. This term brings 
about the nonminimal interaction between the cold DM and chameleon field. 

For a spatially flat FLRW metric variation of the action (I2.90P w.r.t. the 
metric tensor components leads to the following Friedmann equations 


3 H 2 — p m f + -p 2 + V(p); (2.91) 

H 2 + 2- — — i</? 2 + V(p). (2.92) 

a 2 


Variation of the action (12.90(1 w.r.t. field p results in the equation of motion 
for the scalar field 

ip + 3 Hp + V' = —pmf'i (2.93) 

where a prime indicates differentiation w.r.t. p. 

From the equations (l2.91H2.93l) one can obtain the modified conservation 
equation 

j t (p m f) + 3 Hp m f = p m f'f t , (2.94) 

and for the matter density one obtains 


Pm 



(2.95) 


It means that in spite of appearance of interaction with the chameleon field 
the matter density evolution remains unaffected. 

Using the expressions (I2.9ip and (I2.92p . we obtain expressions for p e ff 
and p eff 


Peff — Pmf T T U(^) — Pch T Pdei (2.96) 

Peff = ~ V {p) = Pdei (2.97) 


where de and ch refer to the standard DE and to the chameleon energy re¬ 
spectively, and peff = WeffPeff ■ Here we introduce the so-called ’’chameleon 
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matter density” as p c h = p m f • It worth noting that the density of chameleon 
field, as can be seen from (I2.93p . is characterized by value of p and the poten¬ 
tial V(p), while p c h in its turn describes the DM density interacting with the 
chameleon held. The conservation equations for the DE and the chameleon 
held which interact with matter separately take on the form 


Pch T SHp c h Qi 

Pde T 3H(1 T LOde) Pde Q, 


(2.98) 

(2.99) 


where Q is the interaction source and pd e = WdePde- Comparing 02.941) and 
02.98P one finds that 


Introducing function 


we rewrite Q as 
1 


n - 1 • f _ 1 f 

Q — ^PmPJ — ^ Pch j- 


g(t,p,p) = 


f(<p) 


1 


Q 1 9Pch — gdifi ( I W e//) 2V(1 + UJeff)). 


( 2 . 100 ) 


From equation 02.1001) one observes that t = t cross when u> e ff = — 1, then 

1 .o 


Q = 


:9P 


( 2 . 101 ) 


Current value of the function g directly ahects interaction between the com¬ 
ponents, for instance the interaction is absent at g = 0. 

The equation Q2.100p can be rewritten in the form 


Q = \gpch = -^g(p 2 {2 + q)-2V(l + q )). 


( 2 . 102 ) 


If we consider ratio of the chameleon matter density to that of the DE 
r = p c h/Pde, we then obtain 


1 1 3 H 2 

1 + - =-. 

r Pch 

From equations Q2.103P and 02.98P one finds 

f 


(1 + r) 2 


3rH 2 H Q 
1 + r 1 + 3 H 2 ’ 3 H 2 ' 


(2.103) 


(2.104) 
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In addition, using 


we obtain 



r 

(1 + r ) 2 


tH 
1 + r 


(1 + 2 q) 


Q_ 

3 H 2 ' 


Using equation (12. 1Q0[) now we can find expression for q : 


2(—18# 3 r(l + r) + 18# 2 r + ^(1 + r) 2 (V - <p 2 )) 
(1 + r)(72# 3 r - g( 1 + r)(2U - tp 2 )) 


(2.105) 


(2.106) 


Having four unknown functions (a(t), H(t), <p(t), p m (t )) and only three 
independent equations (12.91112.92112.93]) we are forced to introduce an addi¬ 
tional ansatz. In the paper [53] this ansatz was made in the following form 


H = exp 


1 — 7 a 2 
aa 


(2.107) 


where 7 and a are positive constants. Choice of the ansatz is dictated by 
desire to reproduce well the observed expansion of the Universe. The system 
of equations is now closed, but the high degree of non- linearity forces us to 
make the additional simplifying assumption [53] 

V = (5if 2 . (2.108) 


The condition (3 > 1/2 guarantees negative pressure of the chameleon field. 
In presence of the additonal conditions (12.1071) and (12.1081) we must make 
sure that together with the equations (12.91H2.95]) they provide a consistent 
set of solutions. Such solutions were obtained in [53] 


/ = 


a 3 # 3 


PoW - 1/2) L 

H 2 


(2/3 + 1) ( — + — ] — 3 
' aa a ' 


if 2 = 


-2 \ — 

aa a 


+ 3 


(2.109) 


( 2 . 110 ) 


P - 1/2 L 

Having convinced ourself in existence of the consistent set of solutions, we 
can analyze behavior of the DP consistent set of solutions. 


Q = 


d 

dt 


_ 1 = ± e -na) = e -F(a) (_^La) = 1 + 7 ° 2 

dt V da ) aa 


-i; (2.111) 
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where F(a) = (1 — 7 a 2 )/(cra). We can see that if a > 1 + 7, then at present 
(a = 1) expansion of the Universe is accelerated (q < 0). 

A particular feature of the considered model is non-monotonic dependence 
of DP on the redshift (see Fig J2.4lh 

2.7 Phantom fields 

The enormous number of cosmological observations at our disposal shows 
that the parameter w in the EoS for DE lies in a narrow range in the vicinity 
of w = —1. Above, we examined the region —1 < w < —1/3. The lower 
boundary of the region, w = — 1 , corresponds to the cosmological constant, 
while the entire remaining interval can be realized with the aid of scalar fields 
with a canonical Lagrangian. We recall that the upper boundary w = —1/3 
is related to the necessity of providing the observed accelerating expansion 
of the Universe. Is it possible to go beyond the limits of this interval and 
realize w < — 1? This is a complicated question for the energy component, of 
which we know so little. In GR, certain limits are conventionally imposed on 
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possible values of the components of the energy-momentum tensor by energy 
conditions (see Part I, section 2.6). One of the simplest conditions of this 
type is p + p > 0. The physical motivation for this condition is to avoid the 
instability of the vacuum. Concerning the dynamics of the Universe, this 
condition requires the density of any admissible energy component not to 
increase with the expansion of the Universe. The cosmological constant, for 
which p\ = 0, represents a limiting case. Taking into account our ignorance 
concerning the nature of DE, it is reasonable to ask ourselves: can this 
mysterious substance differ from the ’’good” sources of energy known to us 
and violate the condition p + p > 0? Since DE must have a positive energy 
density (which is necessary for the Universe to be flat) and negative pressure 
(to provide the accelerated expansion of the Universe), this violation should 
lead to w < —1. 

Some time ago, such a component, which was called the phantom energy, 
attracted the attention of physicists [5TM56j . 

The action for a phantom field ip minimally coupled to gravity differs 
from the canonical action for a scalar field by the sign of the kinetic term. 
The energy density and pressure of the phantom field 


Pp = T 00 = + V (<p ), Pip = Ta 




and the EoS has the form 


Pp = 0 2 + 2 V(ip) 
Pp 0 2 -2 V{<p)' 


( 2 . 112 ) 


(2.113) 


If v? 2 < 2V(cp), then w v < —1. As an example, we consider the case when the 
Universe contains only nonrelativistic matter (w = 0) and a phantom field 
(w<p < —1). The densities of these components evolve independently: 


pm CK CL 


and 


Pip CX CL 


3 ( 1 +Wtp) 


If matter domination terminates at a time moment t m , the solution for the 
scale factor at t > t rn is expressed as 


a(t) 


®(tm) 


~VJ<p + (1 + Wp) 



2 

3(1 -\-W(p ) 


(2.114) 
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Hence, it immediately follows that if w v < —1, then, at the moment of time 


t br — 


w u 


(1 + Wp) 


the scale factor and a number of cosmological characteristics of the Universe 
(for instance, the scalar curvature and the energy density of the phantom 
field) become infinite. This catastrophe has been termed the Big Rip. The 
Big Rip is preceded by the so-called superaccelerating expansion. A sim¬ 
ple explanatory example helps to clarify the origin of the superacceleration 
regime. Let us consider the differential equation 


— = Ax 2 . (2.115) 

dt v ’ 

If A > 0 , Eqn (12.115ft realizes a positive feedback. The rapid increase of the 
function x(t) leads to the Big Rip (the function blows up) in a finite time. 
Indeed, the general solution of Eqn (12.115(1 has the form 

x(t) = ——--, (2.116) 

where B is the integration constant. At t — —B , the Big Rip occurs. It is 
easy to see that model (12.1151) represents a concrete version of the Friedmann 
equation for w^ < —1. Because p^ oc the first Friedmann equation 

can be represented as 

a = Aa"i (1+ ^ )+1 . (2.117) 

For example, for = —5/3, Eq. (12.1171) coincides precisely with (12.1151) . 
Note that on the contrary to the phantom dominated case (1 + w v < 0) in 
the quintessence dominated Universe (1 + > 0), the cosmic expansion is 

singularity-free with the scale factor 

a{t) = a (t m ) [1 + (1 + w v ) (t - t m ) /t m \ . (2.118) 


Let us consider for instance a phantom scenario in the case of hybrid expan¬ 
sion law |E7j. In this case the anzatz 


fl(f) — CLq 



(2.119) 
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must be slightly modified in order to acquire self consistency. It can be 
achieved by rescaling the time as t —* t s — t, where t s is a some succinctly 
positive reference time. The rescaling leads to 

a(t) = a °(rr£) (2.120) 

and 

H _ _«_ P_ jj _ a ( __ a ( f s ~ to ) 2 _ 1 

t s -t t s -t 0 (t 8 - ty [/3 (t s - t) + a(t s - t 0 ) 2 ] 2 

(2.121) 

It then immediately follows that a < 0 leads to q < 0 (acceleration). More¬ 
over, as H > 0, then this case realizes the superacceleration q < —1 , H — 
—H 2 { 1 + q). Also the scale factor and Hubble parameter diverge as t —* t s 
and thus exposing the Universe to Big Rip. 

It is also worth noting that in this models at late times the DP tends to 
the constant value 

q(t) ->• ~^\l + 3w\. (2.122) 

The cosmological models with the superacceleration evidently lack big chance 
to be realized in our Universe, as it would make impossible the large-scale 
structures to appear unless too fine tuning of the model parameters is per¬ 
formed. The Universe evolution rather realizes a soft violation of the equality 
w = — 1 , but in such a way that does not (and will not) lead to decay of cos¬ 
mological structures or to the Big Rip. We consider models of such type 

below. 


2.8 Crossing the Phantom Divide 

As was mentioned in the previous section, there exists certain observational 
evidence in favor of the fact that the EoS parameter takes the so-called 
fantom value w < — 1 in the present time or in the recent past. The best 
current measurement [60] for the EoS (assumed constant) is 

w = -1.04±g;?g. (2.123) 

Figure 12.51 shows confidence intervals for the EoS parameter obtained with 
the help of different observational data. 
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w w 

Figure 2.5: Likelihood curves for a constant EoS w in a flat universe, using 
Planck CMB data (left panel) and WMAP9 CMB data (right panel). In 
the figure are compare constraints from CMB + BAO data alone (dashed 
black) to those which additionally include SN la data from SNLS3 (la) |6TJ 
(blue), Union2.1 (green), or PS1 (1.2a) [62J (dashed red). All likelihoods are 
marginalized over other cosmological and nuisance parameters [59] . 


Constraints on the DE EoS from recent SN la data complemented with 
distance measurements from the CMB and a compilation of BAO results. 
For the SNLS3 and PS1 SN data sets, the combined SN la + BAO + Planck 
data favor a phantom EoS where w < — 1 at ~ 1.9a confidence (see Fig. 12.5ft . 
in good agreement with the corresponding results reported in the original SN 
la and CMB data set. 

On the other hand, it is commonly known that the DE with phantom 
EoS violates the dominant energy condition which is a sufficient condition of 
the conservation theorem [53], and therefore the vacuum in such models is 
unstable both on classical and quantum levels. 

Actually, in the articles [Ml 65] they are showing that the scalar fields 
with w < — 1 are quantum-mechanically unstable concerning decay of the 
vacuum into gravitons and phantom particles with negative energy. This 
problem could be solved by introduction of the DE with time-depending EoS 
parameter. If the DE could dynamically change its EoS from a phantom like 
one to that with w > —1, then this crossing would stave off the unacceptable 
particle production. As was noted in Ref. [66JEI], quantum effects in locally 
de Sitter space could lead to the effective parameter w < — 1. 

A further fundamental physical question where this transition could play 
a significant role is the cosmological singularity problem. \i w < —1 in 
an expanding Friedmann Universe, then the positive energy density of such 
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phantom matter commonly becomes infinite in finite time, dominating over 
all other forms of matter and, consequently, leads to the late-time singularity 
called the “Big Rip” [ 68 ] . The models where EoS parameter grows during 
the evolution and the DE transits from the regime with w < — 1 to that with 
w > — 1 naturally prevent this late-time singularity. Here it is worthwhile to 
mention that for certain potentials and initial conditions the phantom scalar 
fields can escape this singularity by evolving to a late-time asymptotic which 
is the de Sitter solution with w — — 1. Moreover, it was argued that the 
quantum effects can prevent the developing of the “big-rip” singularity as 
well. In the present subsection we consider a model of flat Universe filled 
with three components: the DE /?<&, = p\ + p 2 with partial energy densities 
Pi and p2, and non-relativistic matter p w 

H 2 = 7^p(pi + P 2 + p w ) (2.124) 

The EoS for the DE reads 

— = y/A + B(l + z ) 6 , (2.125) 

Po 

— = -y/c + D(l + z) 6 . (2.126) 

Po 

Here p\ denotes the density of type I Chaplygin gas and p 2 represents the 
density of type II Chaplygin gas. 

The continuity equation then is 


dpde = 3 (p de + Pde)d lll(l + z). 


From the continuity equation (12.127ft . we arrive at 


IVde 


Pde 

Pde 


1 d 111 Pde 

3 dln(l + z) ’ 


(2.127) 


(2.128) 


which means that in an expanding Universe if pde decreases and then increases 
w.r.t. the redshift, or increases and then decreases, then we conclude that 
EoS of DE crosses phantom divide. Here, 


dpde 

dz 


—3[Bp 1 X (1 + z) 5 + Dp 2 X (1 + z) 5 ]. 


(2.129) 
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Figure 2.6: Wd e (z ) as a function of z. On all of the 5 orbits B = 0.3, D = 0.15 
with different A and C. From the above to the below A = 1.104, C = 
0.048; A = 1.059, C = 0.039; A = 1.01, C = 0.03; A = 0.9106,0 = 
0.015; A = 0.8267, C = 0.006. 



Figure 2.7: DP q as a function of z. On all of the curves orbits B = 0.3, D = 
0.15, = 0.28 with different A and C. For the solid curve A = 5.92, C = 

3.00, the short dashed curved A = 1.63, C = 0.300, while for the long dashed 
curve A = 1.01, C = 0.0300. 
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Therefore one obtains the EoS parameter as 


w de = V — = -1 + (1 + zf[Bp r 1 + Dpi 1 }. (2.130) 

Pde 


The equation w de (z ) 
-1, w de (z) = -1 


— 1 has a solution. It is easy to see that at z 


Jad 2 -b 2 c 

Ztr ~ _1 + y BD(B — D)' 


(2.131) 


If this expression is treated as a small addition to —1 then one can conclude 
that the regime w de ~ —1 is realized on late stages of evolution of Universe. 
We see that if we carefully tune B and D , the EOS of dark energy may cross 
—1. Figure [2761 shows some concrete examples for which Chaplygin crosses 
the phantom divide at about z = 0.2. 

We now turn to analysis of the DP in this model which is the most 
important parameter from the observational point of view. It is known that 
q carries the total effects of cosmic fluids. Here q reads 


da 1 T . A C . 1 

-T7T = — 1-3— + — - 

o 2 Pi P 2 Pi T P 2 T Pm 


(2.132) 


where pi,p 2 are defined in the expression (12.1251) and (12.126jl respectively, 
and p w is the DE with the EoS parameter w : p w = p w0 (1 + z) 3 * 14 ^. 


q(z = 0) 


1 3 \VaTb Vc+d ) 

2 ~~ 2^ATB-VC + D + 1 


(2.133) 


In the infinitely far future, i.e. at z —> —1, we obtain the de Sitter Universe 
with q —> — 1, and the divergence problem is absent. 


2.9 Quintom Models 

Analysis of properties of the DE from recent observations mildly favor models 
with w crossing —1 in the recent past. However, as was shown above, neither 
quintessence nor phantom can realize this transition. But this transition 
can be easily realized in the Universe where role of the DE is played by 
quintessence and phantom fields simultaneously. This kind of DE is termed as 
quintom [SI] . Besides that, this model for certain parameter values describes 
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the results of observations better than most of other models with w > — 1. 
Moreover, the qnintom DE models provide the so-called tracking field regime 
and therefore facilitate solution of the coincidence problem, in distinction to 
the phantom DE which suffers from this problem particularly heavy [82J. 

In the present subsection we give theoretical basis of the quintom DE and 
also describe some models which reflect its characteristic features. 

As we mentioned before, the quintom model was initially suggested in 
the paper [81], where the authors combined the Quintessence and Phantom 
components together. The action for such a DE reads: 


S, 


Quintom 


= j d 4 xy/=g | i V M <pi V^Vi - - V(<p i, <p 2 ) j , 


(2.134) 

where ip\ and (p 2 are Quintessence and Phantom components, respectively. 

The EoS for such compound DE can provide even multiple crossings of 
the phantom divide due to the combined effects of the two components. In 
general, the corresponding effective potential can have arbitrary form, and 
the fields can have a background interaction as well as a direct one via the 
energy exchange. 

In the original paper [8T] the authors consider a simple model with the 
background interaction: 


_ Ac /?2 

^(^1,^2) = Vo{e m p‘ +e - ™P) 


(2.135) 


where m p i is the Planck mass and A is a dimensionless constant. 

As an example of the Universe with the quintom held we cite the model 
considered in paper [83] . 

Let us consider the simple quintom cosmological model in a hat spacetime. 
The action in such a Universe consists of a canonical 95 and a phantom a 
helds: 


S — J d A X\J—g 


if? - + V a (<x) + ^g^d p ad u a + 14(a) + £ t] 


(2.136) 

where we have set 8ttG = 1. Here V p ((p) is the quintessence, V a (a) is the 
phantom held potentials, while the term C m corresponds to the Lagrangian 
of the (dark) matter content of the Universe, considered as dust. Radiation 
plays significant role on early stages of evolution of the Universe, but we 
consider the later epoch when it can be neglected for the sake of simplicity. 
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The Friedmann equations take on the form: 


— 2 [pm + Pip + Po'j , (2.137) 

H = --(^Pm + + P v + Per + Po^j 1 (2.138) 

and the conservation equations for the quintessence and the phantom fields 
are: 


p<p + 3 H(p lfi + Pp) = 0, 

Pa + 3 H (pfj + Pa) = 0, 

(2.139) 

(2.140) 

In these expressions, p v and p v are respectively the pressure and density of 
the quintessence, while p a and p a are the corresponding quantities for the 
phantom field. Their energy densities and pressures are known to equal 

Pv = \ ( p 2 + v M 

(2.141) 

Pv = 

(2.142) 

Pa = ~^<y 2 + V a (a), 

(2.143) 

Pa = -~a 2 -V a (a). 

(2.144) 

Using the above given equations one can obtain the evolution equations for 
the quintessence and phantom energy: 

p + 3 Hp + = o 

op 

(2.145) 

o + ZHo Sl " (< 7 ) =0. 
da 

(2.146) 


Using the exact solution p m = p m o/a 3 for the DM density, where p m0 is its 
present value, we obtain the closed system of equations. 

In the quintom scenario, the DE field appears as a superposition of the 
quintessence p and phantom field a: 


PDE = Pp + Pa 
PDE = P<p + Pa- 
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Therefore, the DE EoS parameter takes on the form 


Pde Pp + Pa 

w =-= —-. 

Pde + Pa 


(2.147) 


Below we follow the paper (83j to consider the quintom evolution in power-law 
potentials: 


V v {(p) = (2.148) 

V a (a) = K a <J- a % (2.149) 


where and K a are constants with units m 4+QV and m 4+ " a respectively. As 
was shown in the paper (83] the weak energy condition requires a^ to be 
positive and a„ to be negative and bounded, that is potential (12.148ft will be 
an inverse power-law one while (12.149ft will be of a normal power-law form. 

ft is convenient to rewrite the system (12. 145ft - (I2.146ft in terms of the scale 
factoe 


Pp T Pip 

_ 5^ _ !, 

( n dtp \ 2 

(2.150) 

2( Pm. T Pip) 

3 H 2 6 

V G da ) ’ 

Pa+Pa 

~l &2 

- l ( a da '\ 

(2.151) 

2( pm Per) 

3 H 2 

6 V da) 


and thus (l2.141IUI2.143p can be written as: 

\v‘ = TZ~\p™ + VM], (2.152) 

Z 1 X(p 

-\^=T^~W + V a (o)]. (2.153) 

Z i x a 

Therefore, we can simply write: 


Pp> 


Pa 


Pip 


Pa 


•^ipPm T Vipi^p) 

1 - X<p 

£aPm T fA(cr) 
l- x a 

£ipPm 2 x^,) 

1 - X, p 

%aPm f // tr(c r )(l 2 X a ) 
1 ~Xa 


(2.154) 

(2.155) 

(2.156) 

(2.157) 
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and thus for the DE EoS parameter we obtain: 


XipPm. ^X(p ) _i %apm Va-(o")(l ^Xfj ) 

1— X v ' 1 — Xa 

W = --- 

Xippm. + Vpjip) i Xg-Pm + Va-jcr) 

1 — Xp ‘ l—x„ 

Similarly, the first Friedman equation can be written as 


(2.158) 


3 H 2 


Pm + U/(<p) + ^ / o-( Cr ) 
1 - x v - X a 


(2.159) 


Finally, using expressions (12.154^ - 112. 157jl and (I2.159P the field evolution equa¬ 
tions (I2.145P and (I2.146P become 


(2.160) 


(2.161) 

with prime denoting the derivative w.r.t. a. These equations are exact and 
account for the complete dynamics of the quintom scenario. 

Unfortunately there is no exact solution for the system (12.160112.16ip . 
however it is possible to obtain some exact results (see [83] for details) in the 
limit of non-relativistic matter domination, i.e. p^,, \p^\ <C p m and p a , \p a \ <C 
p m , which is relevant only for early Universe, as at the present time the DE 
density is known to be twice greater than that of the DM. We give the results 
of numerical integration of the system (12.160112.16ip . and according to the 
context of our review, instead of the EoS (as it was in [83]) we make accent 
on the DP q, taking into account that the relative density has manifestly 
non-monotonic dependence (see Figure [22]). This is non-trivial distinction, 
at least in the range 1 4 a > 1, for the chosen values of the model parameters. 

DE domination on early stages of evolution of the Universe is charac¬ 
teristic for all DE models. It is the case also for the considered model: in 
the early Universe we get q(a — > 0) —>■ 1/2. While the quintom contribu¬ 
tion in the energy balance grows, this simple picture gets destructed and the 
dependence q(a) acquires clearly non-monotonic character. The first local 


a 2 ip" + ^ (5 - 3x v - 3x a ) + 

3(1 -x v - x a ) f ap'[V v {p) + 14(g)] dV^ip) \ _ 

Pm + V<p(<p) + 14(a) l 2 dp J 

a 2 a" + ( ^~ (5 - ?yx^ - 3x a ) + 

3(1 - Xp - x a ) J- aa'\y^(ip) + 14(a)] _ dV a (cr ) ) _ 

Pm + Vp(<p) + 14(a) I 2 da I 
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Figure 2.8: Relative densities Q m , fide and the DP q as a function of a. All 
curves correspond to k\ =2, k 2 = 0.5, a = —2, f3 = 3 p m = 1 

minimum appears at a ~ 0.862, q fa 0.0135, then follows a short phase of de¬ 
celerated expansion, and the next local minimum at a fa 1.39. q fa 0.2. After 
several pulsations of relative densities (and respective oscillations of the DP 
with general trend to accelerated expansion), relative contribution of the DE 
starts to dominate, and the Universe enters to the stage of accelerated (and 
soon the super-accelerated) expansion. It worth noting that this model is 
free of the Big Rip in spite of the fact that during a long period of evolution 
of the Universe the equivalent EoS parameter lies within the phantom zone, 
and the Universe is in the super-acceleration regime. The minimum value of 
the DP for the given values of the model parameters equals to q m \ n ~ —1.35, 
after which it starts asymptotic approach to the value q(a —> oo) —> —1. 

2.10 Interaction in the dark sector 

SCM considers DM and DE as independent components of energy budget 
of the Universe. DE in the SCM is postulated in form of the cosmological 
constant, introduced as early as by Albert Einstein in order to make possible 
the creation of stationary Universe model. The assumed absence of interac¬ 
tion between the two components means that the energy densities for each 
component obey independent conservation equations 

Pi + 3 H {, pi + = 0 . 

Coupling between the two components leads to modified evolution of the 
Universe. In particular, the energy density for non-relativistic component 
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will not evolve according the law a -3 , and the DE density (in form of the 
cosmological constant) will not remain constant any more. From one hand, 
such modification of the theory opens for us new possibilities for solution of 
principal problems in cosmology. So, for example, solution of the coincidence 
problem!)] reduces to appropriate choice of the interaction parameters 
which can satisfy the condition 


a 


de 


n 


dm 


= 0 ( 1 ) 


in the present time consistently with the condition a > 0. From the other 
hand, introduction of the interaction will modify the relations between the 
observable parameters. In particular, the modification will affect the fun¬ 
damental relation between the photometric distance and the redshift, which 
the evidence of accelerated expansion of the Universe is mainly based on. It 
imposes strict limitations both on the form of the interaction and its param¬ 
eters. 

Taking into account that the DM and the DE represent dominant com¬ 
ponents in the Universe, from point of view of the held theory it is naturally 
to consider interaction between them [861 - 191] . Non-minimal coupling in 

the dark sector can considerably affect history of the cosmological expansion 
of the Universe and evolution of the density fluctuations, thus modifying the 
rate of the cosmological structure growth. Different and independent data of 
many observations, such as Wilkinson Microwave Anisotropy Probe, SNe la 
and BAO, were specially analyzed in order to investigate the limitations on 
the intensity and form of the interaction in the dark sector. Some researchers 
also suggest pm, that dynamical equilibrium of the collapsing structures, 
such as clusters, will essentially depend on form and sign of the interaction 
between the DM and the DE. 


2.10.1 Model of Universe with time dependent cosmo¬ 
logical constant 

The simplest example of the model with interacting DM and DE is the cos¬ 
mological model with decaying vacuum. Actually the A(t)CDM cosmology 
represents one of the cases where the EoS parameter w for DE equals to —1. 


The problem is that during all the Universe history the two densities decay by different 
laws, so it is necessary to impose very strict limitations on their values in early Universe 
in order to make them to be of comparable order nowadays 
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This model is based on the assumption that the DE is nothing that physi¬ 
cal vacuum, and energy density of the latter is calculated on the curved space 
background with subtraction of the renormalized energy density of physical 
vacuum in the flat space [97]. The resulting effective energy density of phys¬ 
ical vacuum depends on space-time curvature and decays from high initial 
values in the early Universe (at strong curvature) to almost zero magnitude 
in the present time. 

Due to the Bianchi identity, the decay of vacuum must be accompanied 
by creation or mass increase of the DM particles, which is common property 
of the decaying vacuum models, or in more general case, for the models with 
interacting DM and DE. 

Now let us give the formulation of the A(t)CDM model. In the present 
section we will use the system of units where the reduced Planck mass equals 
to unity: Mpi = (87 rG ) -1//2 = 1 . 

For the case of flat Universe described by FLRW-metric, the first Fried¬ 
mann and the conservation equation can be presented in the form 

Ptot 3 H , Ptot T 8 II (pi 0 t T Ptot) 0, 

where p to t = p m + Pa, and taking into account that p m = 0 and p\ = — p A , 
the conservation equation takes on the form: 

pm + 3 Hp m = -p A . (2.162) 

The right hand side of the latter contains an additional term which plays role 
of a source — it is the decaying cosmological constant. The problem with 
such approach is in fact that we end up with the same number of equations, 
but acquire additional unknown function A(t). The above described approach 
to define the physical vacuum density, though intuitively clear, faces certain 
principal difficulties [83]. That is why a phenomenological approach prevails 
in the literature. Below we present a simple but exactly solvable model. 
Consider the case when A depends on time as [98] : 

A = aH. (2.163) 

It is interesting to note that with the choice cr m ^ (m n is the energy scale 
of the QCD vacuum condensation) the relation (12.1631) gives the value of A 
which is very close to the observed one. 

For the considered components the first Friedmann equation takes on the 
form 

p 7 + A = 3H 2 . (2.164) 
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The equations (I2.162p and (I2.164p together with the conservation equation 

Pi = w Pl = (7 - 1 ) P 7 > 


and the decay law (I2.162P for the cosmological constant completely determine 
time evolution of the scale factor. Combining those equations, one obtains 
the evolution equation in the following form 

2 H + 3 7-ff 2 — a'yH = 0. 


Under the condition H > 0 the solution for the latter equation reads 

a(t ) = Cfexp ( a'yt/2 ) — l] 3 ^, 


where C is one of the two integration constants ( [109] for details). 
This expression can be used for calculation of the DP 


q(z) = 


f^mo(l + z) 3 ^ 2 


1 — + O m o(l + Z) ! 

One thus finds for the current value of the DP 


- 1. 


(2.165) 


q(z = 0) = -D m0 - 1. 


Therefore the accelerated expansion of the Universe occurs under the follow¬ 
ing condition 




mO < ^- 


2 

3’ 


This condition is in fact satisfied for the observed value D m o ~ 0.23. From 
(12.1651) it follows that the transition from the decelerated expansion to the 
accelerated one took place at 


z 


* 


1-0 


m 0 


o 


raO 


2/3 


-1. 


This value ( z* ~ 1.2) exceeds (being however of the same order) the corre¬ 
sponding value for the SCM (z* ~ 0.75), which is the result of the matter 
production in the vacuum decay process. 

Therefore the model based on the time-dependent cosmological constant, 
where vacuum density linearly depends on the Hubble parameter, appears 
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to be quite competitive. It sufficiently accurately reproduces the ’’canonic” 
results, relevant both to the radiation dominated and matter dominated eras. 
The present Universe expansion is accelerated according to the model under 
consideration. Verification of the model with the latest observational data 
obtained for SNla leads to the results (e.g., 0.27 < f2 m0 < 0.37), that agree 
very well with the currently accepted estimates for the parameters of the 
accelerated expansion of the Universe. 

2.10.2 Varying A as a function of the scale factor. 

Let us now consider a flat two-component Universe filled by matter with the 
EoS p = wp (or p = (7 —1 )p in terms of 7 = tc + 1) and cosmological constant 
with the following scale factor dependence 

A = Ba~ m . (2.166) 

The conservation equation and the first Friedman equation can be used to 
obtain the ODE 

sA“ 37 ) = (£I)”™- < 2 - 167 > 

Its integration leads to the energy density of matter as a function of the scale 
factor 

p(a) = p 0 a~ 3l f(a ), (2.168) 

where as usual a(t 0 ) = a 0 = 1, p(a 0 ) = p 0 , and the function /(a) takes on the 
form 

f - 1) . f „ , 0 

/(o) = l + KoX j 37 -m lf rn A 87 , (2.169) 

{ 37 ln(a) if m = 37 , 

where Ko = B/87iGpo. If rn — 0 then /(a) = 1 and the equation (12.1681) 
recovers the usual result p ~ a~ 3 for the non-relativistic matter (7 = 1) and 
a -4 for the radiation (7 = 4 / 3 ). The new parameter Kq can be determined by 
the decay law (12.16611 and astronomical observations imply B = A 0 = 377q Ao, 
then kq = Ao/Do- 

Substituting (12.1 66 jl and (I2.168jl into the Friedman equation one obtains 
^ = a [Q 0 a~ 3 y(a) - (fi 0 + A 0 - l)a" 2 + A 0 a~ m ] 1/2 , (2.170) 
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where the function /(a) is determined in (12.1691) . The second Friedman 
equation can be rewritten in the form 


^ = (i - y) fi 0 a 1 - 37 /(a) + Aon 1 "™ (2.171) 


Substituting this derivatives into the definition of the DP, one obtains 


da 



(1 - f) npa^fia) + \ 0 a 1 ~ m 
a [O 0 a~ 3 ^f(a) — (D 0 + A 0 — l)a -2 + A 0 a _m ] 


(2.172) 


Depending on values of the model parameters it is possible to obtain qualita¬ 
tively different types of evolution of the Universe (see Figure I2.10.2p . In the 
case 77i > 1, the Universe enters into the decelerated expansion stage starting 
from the z = 0, and in the infinitely distant future expansion of the Universe 
stops 

q(m > 1, z —> —1) —* oo. 

For m = 1, we obtain accelerated expansion of the Universe in the infinitely 
distant future 

. . 1 

qyrn = l,z -> -1) -> 

A distinguishing feature is the finite value of the DP at the distance future: 
for m — 1; 3 we have q(z —> —1) —>■ —1, for m = 3 we have q(z —» —1) —> 5 
and q(z —> —1) — > 1 for m — 2. We note that, the appropriate choice of the 



Figure 2.9: The q(z) dependence for the model (I2.166j) . On the left: for 
D 0 = l 5 7 = l,A 0 = l/2, and m = 1 (dashed), m = 2 (solid), m = 3 (dash- 
dot). On the right: for 7 = 4/3 and m — 1 (solid); m = 2 (dash-dot), m = 3 
(dashed). 
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parameter values, essentially of the parameter Ao which is responsible for the 
decay intensity, would make the DP to take only finite values (at least until 
m < 3). 

Thus, in this simple model there is a large freedom for the fitting. As 
might be expected, with the expansion of the Universe the contribution of the 
cosmological constant is reduced and the rate of the accelerated expansion 
of the Universe also changes, which is reflected by changes of the DP. There 
are three different types of the distant future: eternal slowing down, eternal 
speedup or the so-called stopping of the Universe, i.e. expansion of the 
Universe with ever increasing value of the DP: 

q(z —1) —* oo. 

2.10.3 Interacting DM and DE 

One of the most interesting properties of the DM is its possible interaction 
with the DE. Although the most realistic models (the SCM in particular) 
postulate that the DM and DE are uncoupled, there is no serious evidence 
to consider such assumption as a principle. In the present time many active 
researches are provided to investigate the possible consequences of such inter¬ 
action [MllSGl [941199] 1103] , As we have mentioned above, the interaction of the 
dark components can at least soften some sharp cosmological problems, such 
as the coincidence problem for example. The DE density is approximately 
three times as high as the DM one. Such coincidence could be explained if 
DM was somehow sensitive to DE. 

Remark that the possibility of interaction between the DE in form of 
scalar held and the DM lies in the basis of the warm inflation model. Unlike 
the cold inflation scenario, the former does not assume that the scalar held 
in the inflation period is isolated (uncoupled) from other helds. That is why 
instead of the overcooled Universe in the inflation period, certain quantity 
of radiation is always conserved, which is sufficient to manifest itself in the 
post-inflation dynamics. 

Interaction of the scalar held quantum, responsible for the inflation start 
- the inhaton - with other helds imply that its master equation contains 
the terms describing the process of energy dissipation from inhaton system 
to other particles. Barera and Fang |104j initially assumed that consistent 
description of the inhaton held with energy dissipation requires the master 
equation in form of the Langeven equation, where the fluctuation-dissipation 


56 






relation is present, which uniquely links the held fluctuations and the dis¬ 
sipation of energy. Such equation lies in the the basis of description of the 
warm inflation process. 

Interaction between the components in the Universe must be introduced 
in such a way that preserves the covariance of the energy-momentum tensor 
TdU = 0, therefore T '‘h = — T . v + 0, where u v is the 4-velocity. The 
conservation equations in that case take on the form: 


u T /J,y = —u T ^ = —O 
au± DM\V U v 1 n K :u v- 


de 


For the FLRW-metric the equations (I2.173P transform to: 


Pm "F 3 Hp m Q : 

Pde T 3hf/3 D£; (l + w DE ) Q, 


(2.173) 


(2.174) 

(2.175) 


where p m and p DE are densities of the DM and the DE respectively, w DE is 
the EoS parameter for the DE, H = a/a is the Hubble parameter. If Q < 0, 
then DM continuously decays into DE, and vice versa if Q > 0. Note that 
the equations (12. 174)1 and (12.1751) can be combined to obtain the conservation 
equation: 

Ptot + 3 H ptot(ptot + Ptot) — 0, (2.176) 

where p tot = p DE + p m is the total energy density. 

The interaction between the DM and the DE is effectively equivalent to 
modification of the EoS for the interacting components. Indeed, the equa¬ 
tions (12.1741) can be transformed to the standard form of the conservation 
law written for uncoupled components: 

Pde "F 3i7/9 DB (l + 0 

PdM "F 3 HP DM ,eff O’ 


where 


Q 

W DE,eff W DE 3 ]jp 1 W DM,eff 


Q 

3Hp DM 


(2.177) 


play role of effective EoS for the DE and the DM respectively. 

As we do not know the nature of both the DE and the DM, we cannot 
derive the coupling Q from the first principles [105] . However it is clear from 
the dimension analysis that this quantity must depend on the energy density 
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of one of the dark component, or a combination of both components with 
dimension of energy density, times the quantity of inverse time dimension. 
For the latter it is natural to take the Hubble parameter H. 

The following forms of Q are the most often used in the literature: 


Qi = SHyp m , Q u = 3H-fp DB , Qjjj = 3H'y(p m + p DE ). (2.178) 

We consider below the simplest case Qj\ 

pm T 3 Hp m 5Hp m , 

Pde "k 3 H {p DE T Voe) &Hp m , 


where 5 is the dimensionless coupling constant. 
Integration of the latter equations yields 


Pde 
P m : 


■ p a . 

Pde o u i 

-SPpEu a -[3(l+w nE )+i5] , 
3w de+^ 


^PmO T 


j PPEO A 
3 w de +S J 


3 


(2.179) 


Inserting the obtained densities into the first Friedmann equation and trans¬ 
forming from scale factor to redshift, one gets for H 2 (z ) the following 

n 2 = ir\ + ~ >3 fi + ji m (3«, DE + *)]. (2.180) 

3(3 w DE + d) 


Note that in the considered case the DP q(z) can be easily determined by 
the following 


q(z ) 


1 + zdH 2 
2 H 2 dz 


and direct substitution of (12.1801) results in the explicit dependence of the 
DP on the redshift: 


_ 3 w de Q de (3(1 + w DE ) + 5)(1 + zY 3w de+ 5 ) + Q m + S/w de 

2 3w de Q de (1 + zY 3w oe+ s ) + n m + S/w DE 

(2.181) 

With w DE = — 1 and 5 = 0, this expression coincides with the DP value 
obtained in the SCM. 

Let us now consider the inverse problem: instead of the interaction we 
postulate the ratio 

— = /(a), (2.182) 

Pde 
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where /(a) is arbitrary differentiable function of the scale factor. The inter¬ 
action then takes on the form (see [109] for details): 


fa 


Q E p DE El m I 3w DE ] Hp m Q DE I 3w 


f 


fa 


f 


(2.183) 


Let us define /(a) as 

/(a) = / 0 a«, (2.184) 

where £ is a constant and /o is the present value of /(a). The integration 
constant is png] 


Ptot — Ptotp a 3 [n m o + (1 — f2 m o) a£ ] . 

Inserting this expression into the Friedmann equation, one finds: 

E 2 = ^ = a~ 3 [ Ono + (1 - Ono) ] ~ 3Wxn 
= (1 + z f [^mO + (1 ~~ ^mo) (1 + z) ■ 


(2.185) 


(2.186) 


Using the above obtained formulae, it is easy to find the following expression 
for the DP in the considered model 


Q 1 T 2 de^de ~k w m^m) i 
where relative densities VL de , VL dm have the following form 


n ° E ~ a 


a 


n 


u DEO "L (1 ^DEo)a^ 

and we finally obtain 

1 3 w DE n DE0 (l + zf 


(! ~ to D Eo) a * 

30 + (1 - ^ DE o)a^' 


q = 


2 1 - ^DEO + ^D£o(l + Z f’ 


(2.187) 


(2.188) 


(2.189) 


In the considered model, unlike the SCM, the DE density is not constant, 
so it facilitates solution of the coincidence problem under condition £ < 3. 
The accelerated expansion phase, defined by the relation q(z*) = 0, starts at 
the redshift 

1 

*- = I ,. x 1 7 i ko )' - 1- (2.190) 

(1 + 3w de )Q> deo J 
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This relation naturally generalizes the corresponding result for the SCM: 


V ' ‘mO / 

The difference between the values z* in (I2.190p and (I2.19ip depend on how 
much the parameters £ and w DE differ from 3 and —1 respectively. In the 
case of zero difference the quantities (12 .191 p and (12. 190p coincide. Note that 
at z —y oo Universe expands with deceleration 

1 

9 “> 2 ’ 

the case z —» — 1 corresponds to 

1 2 

q ^2~3 WDE - 

Therefore, as expected, dynamics of the considered model is asymptotically 
(at z —* oo and z —>• —1) identical to the case of the two-component Universe, 
filled by the non-relativistic matter and the DE with the EoS p = w DE p. 
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